N ASA/TM-2003-21 1185 



On the Propagation of Plane Acoustic Waves in a Duct 
With Flexible and Impedance Walls 


Abdelkader Frendi 

University of Alabama in Huntsville 

Huntsville, Alabama 

Bruce Vu 
NASA YA-C2-T 

John F. Kennedy Space Center, Florida 


March 2003 



The NASA STI Program Office. . . in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space sci- 
ence. The NASA Scientific and Technical In- 
formation (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA’s scientific and technical information. 
The NASA STI Program Office provides access 
to the NASA STI Database, the largest collec- 
tion of aeronautical and space science STI in the 
world. The Program Office is also NASA’s in- 
stitutional mechanism for disseminating the re- 
sults of its research and development activities. 
These results are published by NASA in the 
NASA STI Report Series, which includes the 
following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compila- 
tions of significant scientific and technical 
data and information deemed to be of con- 
tinuing reference value. NASA’s counter- 
part of peer-reviewed formal professional 
papers but has less stringent limitations on 
manuscript length and extent of graphic 
presentations. 

• TECHNICAL MEMORANDUM. Scien- 
tific and technical findings that are prelimi- 
nary or of specialized interest, e.g., quick re- 
lease reports, working papers, and bibliog- 
raphies that contain minimal annotation. 
Does not contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored con- 
tractors and grantees. 


• CONFERENCE PUBLICATION. Col- 
lected papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, of- 
ten concerned with subjects having substan- 
tial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA’s 
mission. 

Specialized services that complement the STI 
Program Office’s diverse offerings include cre- 
ating custom thesauri, building customized data 
bases, organizing and publishing research re- 
sults . . . even providing videos. 

For more information about the NASA STI Pro- 
gram Office, see the following: 

• Access the NASA STI Program Home Page 
at http:www.sti.nasa.gov 

• E-Mail your question via the Internet to 
help @ sti.nasa.gov 

• Fax your question to the NASA STI Help 
Desk at (301) 621-0134 

• Telephone the NASA STI Help Desk at 
(301) 621-0390 

• Write to: 

NASA STI Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076-1320 



NASA/TM-2003-21 1185 



On the Propagation of Plane Acoustic Waves in a Duct 
With Flexible and Impedance Walls 


Abdelkader Frendi 

University of Alabama in Huntsville 

Huntsville, Alabama 

Bruce Vu 
NASA YA-C2-T 

John F. Kennedy Space Center, Florida 


National Aeronautics and Space Administration 

John F. Kennedy Space Center, Kennedy Space Center, Florida 32899-0001 


March 2003 



ABSTRACT 


This Technical Memorandum (TM) discusses the harmonic and random plane acoustic waves 
propagating from inside a duct to its surroundings. Various duct surfaces are considered, such as 
rigid, flexible, and impedance. In addition, the effects of a mean flow are studied when the duct 
alone is considered. Results show a significant reduction in overall sound pressure levels 
downstream of the impedance wall for both mean flow and no mean flow cases and for a narrow 
duct. When a wider duct is used, the overall sound pressure level (OSPL) reduction downstream 
of the impedance wall is much smaller. In the far field, the directivity is such that the overall 
sound pressure level is reduced by about 5 decibels (dB) on the side of the impedance wall. 
When a flexible surface is used, the far field directivity becomes asymmetric with an increase in 
the OSPL on the side of the flexible surface of about 7 dB. 
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1. INTRODUCTION 


Low noise emission has become the main focus of all modem machinery. For launch vehicles, 
the problem is even more severe as the vibroacoustic environment is very harsh on all 
surrounding structures. Current launch platform technology is very cumbersome and expensive to 
maintain; therefore, there is an urgent need for simpler and less expensive platforms with faster 
turnaround. In order to achieve this goal, a better understanding of the vibroacoustic environment 
during launch is needed. To this end, the Launch System Testbed at NASA Kennedy Space 
Center has designed a series of experiments 1 to test new launch pad technologies. 

One of the ideas being considered is to direct the exhaust plume into ducts in order to shield the 
launch vehicle from the initial blast wave and the harsh vibroacoustic environment. This will 
eliminate the expensive and cumbersome water injection technology currently in use. In this TM, 
an attempt is made at understanding the propagation of plane acoustic waves over various duct 
surfaces in the presence and absence of a mean flow. Numerical solution of acoustic problems 
requires a high order of accuracy of the numerical technique being used. In this study, the 
modified MacCormack scheme (Gottlieb and Turkel 2 ) is used. The scheme is second-order 
accurate in time and fourth-order accurate in space. This level of accuracy is known to be 
sufficient for the resolution of acoustic problems 3 . In order to resolve sharp gradients without 
spurious oscillations, the fourth-order dissipation of Jameson et al. 4 is used. 

In addition to the accuracy of the scheme, the correct and nonreflective implementation of the 
boundary conditions is critical. In this case, all boundary conditions are derived using the method 
of characteristics 5,6,7 and the far field radiation boundary condition of Bayliss and Turkel. 8 
Ozyoruk and Long 9 showed that a frequency-dependent impedance boundary condition could be 
used in the time domain. Their results were in excellent agreement with impedance tube 
measurements 10 . In this TM, a more generic nonfrequency-dependent impedance boundary 
condition is used. When a flexible surface is used, a fully coupled approach is used to couple the 
vibration of the surface to the surrounding fluid n ’ 12,13 . Full coupling is very important in this 
case, as the pressure loading on the surface is very high. 

2. ACOUSTIC MODELS 

2.1 FORMULATION 


The problem being addressed is that of plane acoustic waves propagating in a duct with various 
wall boundary conditions; i.e., rigid, impedance, and flexible walls. Two computational domains 
are used - one representing the duct alone and the other the duct and its surrounding medium. In 
the case of the duct alone, the mean flow effects are included as in 9 . The governing equations for 
the fluid medium are given by the nonlinear Euler equations written as 


3£ + 3£: + 3£ = o 

dt dx By 


( 1 ) 


where Q is the vector ( p,pu,pv,e) T with p being the density, pu and pv the x and y 
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momentum fluxes, respectively; and e the total energy per unit volume given by 


e = 0.5p(u 2 +v 2 )+ pc v T . 
InEq. (1), the functions F and G are; 
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In addition to Eq. (1), an ideal gas state equation is used 

P = PRT (4) 

where p is the pressure, R the gas constant, and T the temperature. 

When a flexible surface is used on one of the duct walls, the motion of the surface is governed by 




(5) 


where t] is the plate transverse displacement, p p the plate material density, h the plate 

thickness, and T the structural damping. In Eq. (5), D = Eh 2 /12(1- v 2 ) is the plate stiffness 
withE being the modulus of elasticity and v the Poisson ratio of the plate material. The 
coefficient N x in Eq. (5) is given by 


N, 


Eh 
2 L 
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which represents the tension created by stretching of the plate due to bending. In Eq. (6) x 0 is the 
origin of the flexible plate and L its length. The forcing term on the right-hand side of Eq. (5) is 


Ap = p^' — p 


_ outside 


where p is the pressure along the interior wall of the duct and p 
ambient surrounding fluid, respectively. 


(7) 
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Equations (l)-(7) are written in a nondimensional form using the following reference quantities 
for the different variables; 
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The free stream properties used are those of air at ambient conditions; temperature = 289 K, 
density p„ =1.23 kg/m 3 , pressure = 1.013 x 10 s N/m 2 , and sound speed c„ =340 m/sec. 

The ratio of specific heats is y = c p /c v =1.4 and the reference length is l ref = 0.305 m. 


2.2 METHOD OF SOLUTION 

The unsteady Euler equations (Eq. (1)) are solved using an explicit finite difference scheme. The 
scheme, which is a generalization of MacCormack’s scheme obtained by Gottlieb and Turkel 2 , is 
fourth-order accurate in space and second-order accurate in time. The numerical scheme, applied 
to a one-dimensional equation of the form 


du = dF 
dt dx ’ 


(9) 


consists of a predictor step given by 


* 
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followed by a corrector step of the form 


( 10 ) 
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In the above equations, the subscript i denotes the spatial grid point and the superscript n the 
time level. Superscript (* ) corresponds to an intermediate time level (between n and (n + 1) ). 
The fourth-order accuracy is obtained by alternating the scheme given above with its symmetric 
variant. Operator splitting is used to reduce the two-dimensional problem to a sequence of one- 
dimensional problems. If L x and L y denote the solution operators for the one-dimensional x 

and y problems. The solution to Eq. (1) is then obtained by 

(i2) 


3 


The beam equation, given by Eq. (5), is solved using a numerical technique derived by Hoff and 
Pahl 14 . 

In addition to the high-order accuracy required by acoustic problems, the use of nonreflective 
boundary conditions is critical. Any boundary reflection would lead to spurious oscillations and 
nonphysical results. In order to minimize the effect of boundaries on the solutions presented in 
this TM, large computational domains are used whenever possible. In addition, all the boundary 
conditions are specified according to the four characteristics 

p + pcu , p - pcu , v and c„p- p/c„ in the x-direction 

p + pcv, p-pcv, u and c^p-plc^ in the y-direction. (13) 

Over the rigid walls, the normal velocity is zero; whereas, over the flexible surface the normal 
velocity of the fluid is set equal to the out-off-plane velocity of the surface. The four characteris- 
tics of Eq. (13) represent the acoustic waves (first and second), the vorticity wave (the third), and 
the entropy wave (the fourth). In addition to the boundary conditions given by Eq. (13), in the 
far-field, the Bayliss and Turkel 8 radiation boundary condition given by 

P, +c„p r +(p-pj/(2r) = 0 (14) 

is used. This boundary is known to be effective and accurate in simulating outgoing waves. The 
radius r in this case represents the distance from the exit duct axis to the boundary point being 
considered. When simulating an acoustically absorbing surface, i.e., an impedance boundary, the 
normal velocity is no longer zero but is rather computed by the flow solver. For time-harmonic 
waves, the linear, normal momentum equation can be written as 

jiop^v + = 0 . (15) 

Using the definition of impedance given by Z ac = plv leads to 

j(Qp„p + Z ac p y = 0. (16) 

Taking the acoustic impedance of an orifice to be approximately 

vp.-j^T (H) 

AjcxI 

with the denominator representing a surface area (S - Axxl) of the orifice. For a square grid. 
Ax = Ay , the equation satisfied by the pressure at the impedance wall reduces to 

p + 2p y =0. (18) 

The plane acoustics waves introduced at the inflow of the duct are of the form 
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( 19 ) 


P = £ Yj cos(fl> ( f + 0,) 

i 

with £ being the amplitude of the waves and (ct ) l , (j> i ) the frequency and phase. A discrete 
number of frequencies are used (500, 1000, 1500, 2000, 2500, 3000 Hz) with corresponding 
phases of (0, n/3, 2n/3, n , 4n/3, 5nl3). 

When the duct alone is used in the numerical experiments without the outside domain, a mean 
flow of the form 


U(y) = 4M pk c, 



( 20 ) 


is used, with M pk being the peak Mach number and H the duct height. 

3. RESULTS AND DISCUSSIONS 

As previously mentioned, two computational domains were used - one composed of the duct 
alone and the other composed of the duct and its surrounding medium. The two computational 
domains are shown on Figure 1(a) and (b). In the duct alone case (Figure 1(a)), the computational 
domain is 2.75 feet by 0.167 foot in x and y directions with a respective number of points of 
151 and 51. This resolution leads to 15 points per wavelength at the high-frequency of 4000 Hz, 
which is adequate for the high-accuracy scheme used. In the case of duct and its surrounding me- 
dium (Figure 1(b)), the computational domain used is 6 feet by 12 feet in x and y directions 
with a respective number of points of 301 and 601. The duct height is taken to be 1 foot leading 
to 51 equally spaced points in the duct. The per-wavelength resolution is kept the same in both 
domains. 

At first, code validations are performed. To this end, the results obtained by Longatte et al. 15 , for 
a duct with rigid walls, are used. In their numerical experiments on acoustic propagation in two- 
dimensional sheared ducted flows the following inflow disturbance is used 

«(y,r) = £tf(y)sin<*-£)sin(%-f) (21) 

where U (y) is given by Eq. (20), e is the disturbance amplitude and / f { c with / being 

the excitation frequency and = mc„ / 2 H the m lh cutoff frequency of the duct. All the com- 

parisons are made iotM pk = 0.08 and e = 0.05 .For Q = 2, which is below the first cutoff fre- 
quency // , the plane acoustic mode is dominant as shown by Figure 2. Figure 2(a) shows the 
pressure contour plots in the duct and Figure 2(b) shows the time history of the fluctuating pres- 
sure at x = A along the centerline and at the top wall. The pressure contours are indeed those of 
plane waves and the two time histories are identical. 
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For £2 = 8 both plane and oblique waves propagate down the duct as shown by the complex 
pressure contours of Figure 3(a). In addition, the time histories of the fluctuating pressure at 
x- A along the centerline and at the wall are now different as shown by Figure 3(b). For 
£2 = 2 n , the pressure distribution in the duct is shown on Figure 4(a). The figure shows an inter- 
action between a longitudinal and a lateral mode. The time history, Figure 4(b), shows the center- 
line pressure to be higher than that at the wall. All these results are in excellent agreement with 
those of Longatte et al . 15 

3.1 PLANE ACOUSTIC WAVE PROPAGATION IN A DUCT WITH AND WITHOUT 
MEAN FLOW 

3.1.1 M pk =0 

Numerical experiments are carried out to study the effects of the various surfaces on the 
propagation of plane acoustic waves, obtained using Eq. (19), in a duct. The length of the 
impedance wall is 38 centimeters (cm) while that of the flexible surface is 28 cm. The properties 
of the flexible surface are stiffness D = 40 Nm, mass per unit are p p h = 2.714 kg/m 2 m, 

Poisson’s ratio v = 0.3, physical damping T = 1.52xl0 2 Ns/m 3 , and a length L - 0.304 m. The 
first five natural frequencies are 134, 371, 727, 1198, and 1783 Hz. Figure 5 shows the pressure 
contour plots for the various surface combinations; rigid top and bottom, rigid top and impedance 
bottom, flexible top and rigid bottom, and flexible top and impedance bottom. The figure shows 
that for the rigid surface case, the pressure contours are vertical, whereas in all the other cases the 
pressure contours exhibit a more complex shape. In particular, one can observe more pressure 
contours over the flexible surface and curved contours over the impedance wall. 

Figure 6 shows the variation of the OSPL over the top wall of the duct as a function of the 
downstream distance for the various surfaces. It is clear from the figure that the case with 
impedance wall exhibits a drop in OSPL of nearly 30 dB near downstream exit of the duct. A 
negligible increase in OSPL is obtained over the flexible surface. When a flexible surface is used 
on the top wall and an impedance surface is used on the bottom wall, the OSPL near the exit of 
the duct is higher by about 5 dB over that given by the impedance wall alone. Therefore, the 
effect of the flexible surface is to increase the OSPL in the duct. Figure 7 shows the effect of the 
impedance wall on individual frequencies. The figure shows that at all the frequencies, the effect 
is nearly the same. This is expected, as the boundary condition used is frequency independent. 
Therefore, for this narrow duct, the impedance wall suppressed most of the acoustic disturbances 
downstream. However, when the duct width is increased from 0.167 foot to 1.0 foot the effect of 
the impedance wall on the OSPL measured on the top wall is negligible, Figure 8. This indicates 
that using an impedance wall to suppress plane waves propagating in a wide duct will not be as 
effective as in a narrow duct. 

Figure 9 shows a comparison of power spectra near the center of the duct’s top wall. The largest 
decrease in sound pressure level at all frequencies is achieved using the impedance wall alone. 
The decrease is more pronounced at 1000 Hz than it is at 3000 Hz at this location. When a 
flexible surface or a combination of a flexible surface and an impedance wall is used, additional 
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frequencies are obtained that correspond to the structural frequencies and the SPL is slightly 
higher. 

Figure 10 shows the time history of the pressure near the center of the top wall of the duct for the 
four different cases studied. In agreement with the previous results, the impedance wall case 
shows the most reduction in pressure fluctuations; whereas in the case of the flexible surface, an 
increase in fluctuations is obtained. Figure 1 1 shows the instantaneous spatial variation of the 
pressure along the duct center for the four different cases studied. Similar to the time history, the 
case of the impedance wall alone results in a large reduction of the amplitude of the pressure 
fluctuations. The flexible surface increases slightly the amplitude of oscillations. 

3.1.2 M pk =0.5 

When a mean flow was introduced, the results changed as follows. The pressure distribution in 
the duct shown on Figure 12 shows less reduction in amplitude of the pressure over an imped- 
ance wall. Also the curved contours are now confined to the region near the wall. The flexible 
surface contribution is also reduced. The variation of the overall sound pressure level over the 
top wall of the duct with the downstream distance is shown on Figure 13. There is a significant 
reduction of the effect of the impedance wall on the OSPL. The decrease in OSPL over the length 
of the duct is reduced to about 20 dB from the 30 dB for the no-mean flow case. When a flexible 
surface is combined with an impedance wall, the reduction is even smaller (around 15 dB). The 
power spectra at the center of the duct’s top wall are shown on Figure 14. The effect of the mean 
flow is to reduce the damping effect of the impedance wall and the contribution of the flexible 
surface. The time history of the pressure at the center of the duct’s top wall is shown on Figure 
15. Similar to the no-flow case, the impedance wall leads to a reduction in the amplitude of the 
pressure fluctuations. The flexible surface, however, has very little effect. The spatial variation of 
the pressure at the duct center for the four different cases is shown on Figure 16. Similar to th.e 
time history, the amplitude of the pressure fluctuation is reduced with the use of an impedance 
surface. 

3.2 DIRECTIVITY OF ACOUSTIC WAVES PROPAGATING AWAY FROM A DUCT 

In this section only the no mean flow results are presented (M pk =0). The duct in this case is 

28-cm wide and 84-cm long. Figure 17 shows a comparison of pressure contour plots for the 
rigid surface case and the case with an impedance wall at the bottom and for a plane wave fre- 
quency of 500 Hz. There are small differences between the two plots both inside and outside the 
duct. Similarly, a comparison of pressure contour plots between a rigid and impedance-flexible 
surface on the top wall of the duct shows a very clear difference, Figure 18. Radiation from the 
flexible surface is dominating the acoustic field both inside and outside the duct. Figure 19 shows 
the directivity pattern along a semi-circle of radius 2 feet centered at the exit of the duct. The use 
of the impedance surface results in slightly lowering the far field OSPL. The strongest effect is 
obtained when a flexible surface is used on the top wall. In this case the directivity pattern is 
asymmetric and an increase of about 7 dB is obtained at various points in the upper quadrant. 
Combining a flexible surface with an impedance wall slightly reduces the OSPL. The time his- 
tory of the pressure 2 feet downstream of the exit of the duct is shown on Figure 20. Slightly 
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lower fluctuation amplitudes are obtained with the use of the impedance wall, whereas the flexi- 
ble surface increases the level of the fluctuations and the frequency content. The spatial variation 
of the pressure along the centerline of the duct and outside is shown on Figure 21. The figure 
shows the propagation and decay of the acoustic waves. Enhanced decay is obtained with the im- 
pedance wall. 

4. CONCLUDING REMARKS 

Results are obtained for plane acoustic wave propagation in a duct over rigid, impedance, and 
flexible surfaces with and without mean flow. The results showed that using the impedance wall 
leads to a decrease in the overall sound pressure level downstream. However, the presence of a 
flexible surface results in an increase of the OSPL in the duct and outside. Downstream of the 
duct, the directivity pattern is also affected by the various duct surfaces. Using a flexible surface 
results in an asymmetric directivity with an increase in OSPL on the side of the flexible surface. 
The impedance wall also affects the directivity, but only slightly. This later result is found to be 
due to the increased duct width. Therefore the impedance wall is more effective at reducing the 
downstream OSPL in narrow ducts. 
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Figure 1. Computational domains: (a) duct alone, (b) duct and surroundings 
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(a) 




Figure 2. (a) Pressure distribution in the duct for ft = 2 , (b) Time history of the pressure fluctua- 
tions at x-k 
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Figure 3. (a) Pressure distribution in the duct for ft = 8 , (b) Time history of the pressure fluctua- 
tions at x = A 
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Figure 4. (a) Pressure distribution in the duct for Q, = 2n , (b) Time history of the pressure 
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Figure 6. Overall Sound Pressure Level (OSPL) along the top wall of the duct 



Downstream distance, ft 

Figure 7. Sound Pressure Level along the top wall for the various frequencies and for the case of 

the impedance bottom with M pk = 0 
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Figure 8. Sound Pressure Level along the top wall for / = 1000 Hz and for the case of the im- 
pedance bottom with M pk = 0 






Figure 9. Comparison of pressure power spectra at the center of the top wall of the duct obtained 

with the various surfaces 
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Figure 10. Comparison of the time history of the fluctuating pressure at the center of the top wall 

of the duct obtained for the various surfaces 


Wall B.Cs. 



Figure 11. Comparison of the instantaneous spatial variation of the fluctuating pressure 

obtained for the various surfaces 
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Figure 14. Comparison of pressure power spectra at the center of the top wall of the duct ob- 
tained with the various surfaces 

WaMB.Cs. WallB.Cs. 




Figure 15. Comparison of the time history of the fluctuating pressure at the center of the top wall 

of the duct obtained for the various surfaces 
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Figure 16. Comparison of the instantaneous spatial variation of the fluctuating pressure 

obtained for the various surfaces 
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Figure 17. Instantaneous pressure contours in the duct and its surroundings for the rigid and 

impedance surface cases (M pk = 0, f = 500 Hz) 
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Rigid Wall Boundary Conditions 


Flexible Wall Boundary Conditions 
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Figure 18. Instantaneous pressure contours in the duct and its surroundings for the rigid and the 
impedance-flexible surface cases (M pk = 0 , / = 500 Hz) 


90 ° 



Figure 19. Far field SPL directivity obtained with the various surfaces at a distance of 2 feet from 

the exit of the duct (M pk =0) 
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